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Abstract 

Combining v oscillations at T2K or T2HK with V oscillations from p + decay at rest (DAR) 
allows a determination of the leptonic CP-violating phase 5. The degeneracies of this phase 
with #i 3 and #23 are broken and 5 can be reliably distinguished from 180° — 6. We present 
the sensitivity to 5 of T2(H)K together with a p + DAR experiment using Super-K as a near 
detector and Hyper-K at the Tochibora site as a far detector. 
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1 Introduction 


We propose the experiment Tokai ’N Toyama to Kamioka (TNT2K). The 50 kton water 
Cherenkov detector Super-Kamiokande (SK) will detect both v e appearance in a n /t beam 
created at J-PARC in Tokai and also v e appearance in isotropic v jL created by // 4 decay at 
rest (/iDAR) at a high intensity accelerator just south of Toyama. We will show that this 
yields a precise determination of the lcptonic CP-violating phase S. With the addition of 
just one fifth of Hyper-Kamiokande (HK) at the preferred Tochibora site, we find that S can 
be determined more precisely than with T2HK using the full megaton HK but no /iDAR 
and also S can be reliably distinguished from 180° — 6. 

The J-PARC beam creates by colliding a 750 kW, 30 GeV proton beam into a target, 
ejecting n + which are filtered into a tunnel where they decay into a /i + and As v 
oscillations are provided by /iDAR from the Toyama accelerator, an optimal determination 
of 6 arises when the J-PARC beam runs exclusively in v mode. The ry travel along the 
beam, oscillating as they go, and SK detects, largely via charged current quasielastic (CCQE) 
interactions, both v e appearance in the beam and also the disappearance. SK is located 
off of the beam axis, where the beam is relatively monochromatic, centered on the first 
oscillation maximum which at 295 km is at about 600 MeV. The Tochibora HK site is at 
the same baseline and off-axis angle, and so if built will also observe v at the first oscillation 
maximum. 

The /iDAR source collides protons into a target, creating 7r + and 7r _ . The proton beam 
must be of low enough energy, and the target large enough, so that the resulting tt stop in 
the target. The tt~ are absorbed or decay into and /i~ which are absorbed in a sufficiently 
high Z target. On the other hand the 7r + decay at rest, yielding /i + and u^. The /i + in turn 
also stop and decay at rest, producing e + , and v e . The spectra of all of these neutrinos 
are illustrated in Fig. 1. The experiment searches for the conversion of to v e between the 
source and the detector(s). The /i + decay at rest spectrum is known quite precisely. Most 
of the v e will have energies of between 30 MeV and 50 MeV and so will interact with SK 
and HK via inverse /3 decay (IBD), whose cross section is also known quite precisely. IBD 
creates an additional neutron whose capture SK-IV is sometimes able to detect [1] and use 
to reduce backgrounds. It is expected that, despite its lower PMT coverage, HK will have 
the same ability. However, to be conservative, we do not use this in our analysis. 

SK and HK can easily distinguish the low energy /iDAR V e from pulsed, higher energy 
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Figure 1: The neutrino spectrum produced by n + and /t + decay at rest (DAR). The charge 
conjugates are produced by tt~ and /i~ decay at rest, which we assume to be suppressed by 
a factor of 5 x 1CT 4 . 


J-PARC v and so the //DAR and beam v experiments can run simultaneously. We will show 
that the optimal location for the //DAR source is 15 km north of SK, which is 23 km north 
of ffK, in the southern hills of Toyama. Fortunately the optimal distance to SK is roughly 
the same while HK is not in operation. Even if ffK is not constructed, with a combination 
of //DAR and the T2K v beam one can determine 5 with moderate accuracy. The favored 
location for ffK is the Tochibora mine, 8 km south of SK [2, 3]. ffowever, a site in the 
Mozumi mine has also been considered. As this location is very close to SK, if it is chosen 
then the TNT2K experiment will be less precise, suffering from the 6 180° — 6 degeneracy. 

The Toyama accelerator needs to accelerate protons to between 400 MeV and 3 GeV, with 
an optimal performance per MW between 600 MeV and 1.5 GeV. So far the most powerful 
such accelerator is the 2.2 mA, 600 MeV proton cyclotron at the Paul Scherrer Institute. As 
explained in Ref. [4] , an optimization of that design would allow for a 5 mA, 800 MeV proton 
beam, which would already be sufficient for our purposes. The target station in, for example, 
Ref. [5] would be suitable. On the other hand, there are currently efforts around the world 
to run accelerator driven subcritical reactors. For example, the project [6] envisages a 10 
mA, 1.5 GW beam with an intermediate stage which is already 10 mA at 800 MeV. For 
concreteness we will choose intermediate parameters, a 9 mA, 800 MeV beam which could 
be achieved for example with two of the accelerators of Ref. [4], 

In order to determine the flux normalization we also require a small, near detector. 
For example a 20 ton liquid scintillator detector, such as one of the eight which Daya Bay 
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will no longer need or perhaps one from RENO, would be quite sufficient. Using elastic 
scattering, such a detector should be able to determine the flux normalization to within 
about 5% [7] and also to provide a very powerful check of the LSND anomaly [ 8 ] with a 
reach to lower mixings and mass splittings than LSND itself. 

At the Erst oscillation maximum, the uncertainty in sin 2 (20i3)sin 2 (023) is about half as 
large as the maximal contribution of sin(5) to the —y v e oscillation probability P^ e and 
so —y v e oscillations alone cannot demonstrate leptonic CP-violating beyond the 2a level. 
This problem can be resolved by combining —» v e and v M —» v e oscillations because, at the 
first oscillation maximum, the sum of the oscillation probabilities P^ e + Ppg depends upon 
sin 2 ( 26 l i 3 )sin 2 ( 023 ) and to a lesser extent on cos(5) while the difference P^ e — Ppg depends 
upon sin(A). Thus by comparing the v beam and /iDAR V experiments one can accurately 
extract sin(5), which measures leptonic CP-violation. 

On the other hand, the existence of multiple baselines is useful not only to control sys¬ 
tematic errors [9] but also to extract cos(<5) [10], thus breaking the S —> 180° — 5 degeneracy 
present in MINOS, T2K, NOz/A, MOMENT and other beam experiments [11, 12], Motiva¬ 
tion for measuring cosfA) is given in Refs. [13, 14, 15, 16, 17, 18, 19, 20]. 

While our proposal for the measurement of 5 is in spirit similar to that of the DAE5ALUS 
project [21, 9], it differs in one key respect. DAE5ALUS uses 3 cyclotron complexes as /rDAR 
sources. However, as the direction of the IBD positron is only very weakly correlated with 
that of the incoming F e , the spectra from the different sources can only be separated by 
running just one accelerator at a time. Thus each runs with a duty factor of only 20% 
and so requires an extremely high instantaneous intensity. To achieve this high intensity 
it will accelerate 77,/" molecules but this involves technological progress, for example, the 
excited molecules must be removed. Our proposal is designed to be cheaper because only 
a single cyclotron complex is necessary and, as it may in principle run continuously, the 
instantaneous intensity may be reduced by up to a factor of five. 

The 20% duty factor at DAEDALUS serves not only so that only one cyclotron runs at 
a time, but also the 40% dead time allows one to measure backgrounds. However, as will 
be explained below, by far the dominant background at our /uDAR experiment arises from 
invisible muons created by atmospheric neutrinos. These lead to a background with a known 
shape, so only the normalization must be determined. However SK has been measuring this 
background, as part of its diffuse supernova neutrino search, for nearly 20 years [22, 23, 24], 
Thus, while some accelerator downtime is inevitable and this will be used to measure the 
atmospheric backgrounds, which are after all dependent on the season and solar activity, we 
do not require a structured beam. 
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In this regard the locations of our experiment in western Japan provide yet another 
advantage. The main backgrounds arise from low energy atmospheric neutrinos. According 
to the model of Ref. [25] at the Kamioka mines the horizontal component of the geomagnetic 
field is 0.31 Gauss, appreciably higher than the 0.17 Gauss that may be expected at DUNE 
or the 0.13 Gauss at LENA in the Pyhasalmi mine. As this strong horizontal field deflects 
low energy cosmic rays, the atmospheric neutrino backgrounds at the sites suggested in our 
proposal will be reduced by nearly one half [26, 27] with respect to the other sites at which 
/i + DAR measurements of S have been proposed. 

2 Parameters 

2.1 The Neutrino Mass Matrix 

We fix the solar neutrino mass splitting to be 

AM 2 2 j = 7.5 ± 0.2 x 10~ 5 eV 2 (2.1) 

and, for ease of comparison with previous studies, we choose the neutrino mass matrix mixing 
angles to be 

sin 2 (20i 3 ) = 0.089, sin 2 (20i 2 ) = 0.857, sin 2 (0 23 ) = ^ (2.2) 

The current uncertainties [28, 29, 30] 

Jsin 2 (20 12 ) = 0.024, Jsin 2 (20 13 ) = 0.005 

Jsin(0 23 ) = 0.055 (2.3) 


are used. 

We assume that the neutrino mass hierarchy has already been determined when our 
experiment has collected its data. While the choice of the hierarchy has little effect on the 
sensitivity to 5, the assumption that it is known does break a degeneracy in experiments, such 
as T2K, NOz/A and T2HK, in which the matter effect is appreciable. As the experiments 
proposed here are unlikely to be performed in the next 10 years, it is reasonable to assume 
that the hierarchy will be known with some cautious certainty. 

Currently, the atmospheric mass splitting A [31] has only been measured at the 4% 
level [32, 30]. Daya Bay has matched this precision [33] for the corresponding effective mass 
splitting AMg e [31]. However, the disappearance channels at T2K and NOz/A will each 
achieve a better than 2% precision, with a 1% precision possible when they are combined 
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[34], while JUNO and RENO-50 are each expected to achieve a subpercent precision [35]. We 
will consider the uncertainty of the MINOS measurement [32] together with the hierarchy- 
averaged central value 

AMjj = (2.4 ± 0.1) x 10~ 3 eV 2 . (2.4) 

Recently SK has reported an excess of —> v e in the J-PARC beam [36] and also a low 
energy atmospheric v e excess corresponding to a deficit in v e —>■ v fl [1], which together give 
roughly a 2<r preference to S — 240° over 6 = 60°, although null CP violation is allowed 
within In. On the other hand, a small deficit in accelerator neutrino —>■ v e at MINOS 
has led to a statistically insignificant preference for sin(J) > 0 [37]. In this note we will not 
consider any of these hints in our analysis. 


2.2 Experimental Setup 


We normalize the detector efficiency and the /i + DAR rate such that, at 10 km if <5 = 0, 350 
inverse /3 decay events, corresponding to v e capture on free protons in SK would be observed 
in a 6 year run. By scaling results from LSND [8] this roughly corresponds to 6 years of 
collisions of an 800 MeV proton beam on a stationary target if the beam current is 9 m A. 
The beam is not pulsed: although clearly a real beam will have dead time which can serve to 
measure the background, we approximate our duty factor to be 100%. Thus the maximum 
instantaneous current is also 9 mA, a factor of 4 less than that which will be required at 
DAEJALUS [38]. The integrated current corresponds to a total of 1.1 x 10 25 protons on 
target (POT), 60 times more than LSND. 

As SK and the Tochibora HK site are only separated by 8 km, the difference in the two 
baselines can at most be 8 km. For the TNT2K experiment we will place the p + source 
15 km north of SK, just south of Toyama city, as is illustrated in Fig. 2. Thus the near 
and far baselines will be 15 and 23 km, respectively. We will assume that J-PARC offers 
750 kW of its beam to this effort, in line with the goal in the next 5 years stated in KEK’s 
most recent road map [39]. Note that this is less than half of the beam power traditionally 
considered in simulations of T2HK [2] . 

We consider SK and HK without gadolinium [40]. As a result of statistical fluctuations 
in the number of photoelectons, in the energy range relevant to the DAR experiments we 
will consider fractional energy resolutions of 


5E 40% , 60% 

—— = — . and — . 

E yB/MeV yE/MeV 


(2.5) 


respectively. 
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Figure 2: Hyper-K and Super-K will be the near and far detectors for p + DAR V from an 
accelerator complex just south of Toyama, while simultaneously detecting higher energy v 
from the J-PARC beam. 


In the case of higher energy accelerator neutrinos from J-PARC, the energy resolutions 
are no longer limited entirely by photoelectron statistics, and so are somewhat worse than 
one would extrapolate from Eq. (2.5) as can be seen in Refs. [41] and [2] for SK and HK 
respectively. In addition, A resonance charged current interactions transfer some of the 
neutrino energy into additional pions and so yield an average energy which is reduced by 
about 360 MeV [42] . We incorporate the reduction in energy of some events and the energy 
resolution by folding the true spectrum with the sum of three Gaussians whose forms are 
given in Appendix A of Ref. [42] . 

The expected p + DAR signals are reported in Fig 3 while the expected appearance 
and disappearance spectra at T2K and T2HK are reported in Fig 4. These were obtained 
using the NuPro package [43] and, in many cases, were confronted with the results of an 
independent C++ code [44], We used definition of T2(H)K in Ref. [45] with a target fiducial 
mass of 22.5 kton (560 kton) and a 750 kW beams operating 10 7 seconds/year for 12 years 
(6 years) in v mode. 

2.3 Flux uncertainty 

The uncertainty in the J-PARC v and V rates (flux times cross-section) is taken to be 5%, 
with the uncertainties uncorrelated. On the other hand, we make a crude approximation 
that they are 100% correlated between T2K and T2HK. 

The uncertainty in the total DAR event rate reflects the uncertainty in the p + DAR rate 
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Figure 3: The expected 12 and 6 year //DAR signals and backgrounds at the 15 km and 23 
km of the TNT2K detectors SK and HK in the case of the normal hierarchy. 


v e at T2K [12 years] 




v e at T2HK [6 years] v^, at T2HK [6 years] 




Figure 4: The expected appearance and disappearance channel spectra at T2K (T2HK) after 
running for 6 (12) years in the v mode. 
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itself at the target and also the efficiencies of the detector. The former can be determined 
using dead reckoning with an accuracy of about 20%. As has been demonstrated by LSND 
[8], a relatively small water or mineral oil based liquid scintillator near detector has a good 
energy resolution in the relevant energy range and can determine the flux much more precisely 
via various channels, such as neutrino electron elastic scattering, which can be separated 
from CCQE interactions using the fact that the angular distribution of elastically scattered 
electrons is strongly forward peaked. We will assume that a near detector is built and so, 
together with a calibration of the DAR, the event rate normalization error can be reduced 
to only 5%. 

One choice for a near detector may be the 50 ton liquid scintillator detector proposed in 
Refs. [46, 47]. The authors proposed that this detector be built 17 meters from a 0.33 mA, 3 
GeV accelerator and search for z7 ;i — y v e oscillations which would indicate sterile neutrinos. 
As this baseline coincides with that of KARMEN, which observed no such oscillations, the 
high AM 2 regime probed by this experiment is largely excluded. The authors mention that 
at a later time it may be desirable to build a farther detector to search the lower AM 2 regime, 
but to obtain a sufficient number of events such a detector would need to be much larger. In 
addition, the 3 GeV beam energy is not optimal for a DAR experiment. It produces less DAR 
events per unit of beam power than a 800 MeV proton beam and has higher backgrounds, 
although the backgrounds can be reduced using the time structure of the 3 GeV beam. 

The present proposal would address both of these issues. The 800 MeV proton beam 
discussed here has 27 times the current of the J-PARC beam of Refs. [46, 47], even after an 
extensive upgrade. Therefore the near detector may be placed 50 meters from the accelerator 
instead of 17 m. This is even longer than the baseline at LSND, and so would extend the 
reach in AM 2 for a sterile neutrino search throughout the entire region suggested by LSND. 
Even at 50 m the detector would see three times as many events as the original proposal 
at 17 m. The improved statistics mean that the shape of the observed spectrum could be 
used to differentiate sterile neutrino oscillations from other potential signals. Furthermore, 
by tripling the distance from the accelerator a number of accelerator-related backgrounds 
are greatly reduced, as well as low energy cosmogenic muon backgrounds if the detector is 
placed 50 meters underground. 

3 Backgrounds 

Our signal and background rates are summarized in Fig. 3 and Table 3. The shape of the 
signal spectrum is well known, it is the spectrum of V M from /x + decay at rest. It vanishes 
above its maximum at 53 MeV and is appreciable above 30 MeV. Below 10 MeV reactor 


neutrino backgrounds hopelessly dominate it. Below 20 MeV it will be dominated by the 
decays of spallation products created by cosmogenic muons, although a double coincidence 
with neutron capture could reduce this background considerably. Near 20 MeV the signal 
may also be subdominant to the diffuse supernova background. The regime from 30 MeV to 
53 MeV will be the main focus of these experiments, containing the first oscillation peak of V e 
appearance. While the second oscillation peak is more sensitive to 6, it would be extremely 
challenging to resolve above the backgrounds. 

Our signal arises from inverse beta capture of oscillated v e , from /r + DAR, on free protons. 
However the target volume also contains oxygen nuclei. These can interact via a quasielastic 
interaction with both the oscillated z/ e and also the much more numerous unoscillated v e 
from /i + DAR. The rate of the CCQE interaction 

v e + 16 0 —> e“ + 16 F (3.1) 

is actually higher than our signal rate. However the Q value of 15.9 MeV, combined with 
the fairly low v e flux above 45 MeV, implies that only a small fraction of the electrons have 
energies above 29 MeV and so will be removed by our low energy veto. This leaves the 
background shown in Fig. 3, which is only appreciable below 35 MeV. It was calculated by 
folding the /iDAR v e spectrum into GENIE and shifting the energy by hand to reproduce 
the correct Q value, and renormalizing the event rate to agree with the calculated values 
in Ref. [48]. Note that this background is much larger at DAEhALUS [21], where the low 
energy veto is 20 MeV. For completeness, in Fig. 3 we have also included the u-e~ elastic 
scattering background, although the vast majority of these events can be removed using an 
angular veto [7]. 

There are also several beam on backgrounds. By far the largest of these arrives as follows. 
When the beam hits the target, not only n + are made, but also tt~. For example, at LSND 
the ratio of n~ to 7r + is about 1 to 8 [7]. Most n~ stop in the target and are immediately 
absorbed. However some, about 5% at LSND, decay in flight yielding n ~. At least 90% of 
fi~ are absorbed in a high Z (Z > 20) target [49]. The rest stop or are captured into orbit 
about oxygen nuclcii. Their decay, which we loosely refer to as /i”DAR, yields an irreducible 
background of V e with a spectrum which is similar to our signal. This background is small 
but difficult to quantify. Roughly following the estimates above, we have fixed the ratio of 
/r - DAR to /i + DAR to be 5 x 10 -4 . This is quite conservative as, at LSND, these decays in 
flight occur in a vacuum region following the target, but we require no such vacuum region. 
In this study we have assumed an uncertainty in this ratio of only 5% (of 5 x 10~ 4 ), however 
in [50] we found that the results are not significantly changed using an uncertainty of 25%. 
We have found that this background is significant only for baselines of less than 5 km. 
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3.1 Atmospheric Electron Neutrinos 

There is an irreducible background due to low energy atmospheric neutrinos. Atmospheric 
electron antineutrinos in the energy range of 30 MeV to 53 MeV IBD capture on hydrogen 
identically to the signal z/ e , although the shape of the background is quite different from that 
of the signal and so a shape analysis can be applied. Atmospheric v e and u e of higher energy 
may enjoy quasiclastic (QE) charged current (CC) interactions with oxygen in the water 
Cherenkov detector. These QE background events often result in the creation of additional 
particles which can be used to veto them. 

We use the unoscillated V and v spectra at the Kamioka mines given in Ref. [27] at 
energies above 100 MeV and at lower energies we use the spectra available on M. Honda’s 
website [51]. Neutrino oscillation is performed for several sample points in each angular bin 
and then angular integration yields the oscillated flux expected at the Kamioka site. 

IBD events are charged current interactions of v e with a free proton which yield a neutron 
and a positron. We will approximate free protons to be 11% of the fiducial mass, and so 
2.4 kton and 62 kton at SK and HK respectively. Therefore there will be 1.4 x 10 33 free 
protons in SK and 3.7 x 10 34 in HK. Multiplying the oscillated neutrino flux, the total 
detector cross section and the 6 -year runtime we find 9 IBD background events between 30 
and 54 MeV at SK and 225 at HK. 

We use GENIE simulations [52] to calculate the number of electron and positron events 
that will result from charged current quasielastic interactions (CCQE) on oxygen. In contrast 
with IBD z7 e , the corresponding v e and V e energies are generally much higher than 50 MeV, 
however due to the Fermi momentum of the nucleon target, the resulting charged lepton 
energy can be in our signal range. We fold the resulting electron and positron spectra with 
the oscillated atmospheric v e and V e spectra to derive the expected atmospheric v e and z7 e 
backgrounds. We consider the sum of the IBD v e and CCQE v e and V e backgrounds. 

SK-IV is able to detect 7 from the H capture of the neutron arising from IBD interactions 
of v [1] with an efficiency which is now about 20%. We find that a double coincidence 
requirement could in principle eliminate most of the CCQE background. However, to be 
conservative we have not applied this veto in our analysis. 

3.2 Invisible Muons 

TNT2K faces an additional, larger, background. Charged current interactions of atmo¬ 
spheric muon neutrinos on oxygen in the detector will produce muons. Those muons below 
the Cherenkov threshold will be invisible to SK and HK and those that decay will produce 
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electrons or positrons whose signal constitutes the invisible muon background. More specif¬ 
ically, Up yield /i~ of which about 20 % will be absorbed by the oxygen in the water and two 
thirds of the relevant events will produce no neutrons and so will not yield a false double 
coincidence. On the other hand produce fi + , which are not absorbed, and usually yield 
neutrons and so a fake double coincidence. While this makes each event more dangerous, 
the higher cross-section for events in this energy range in fact implies that they provide 
the dominant background. While most of these z/ M events can in principle be vetoed by 
requiring a double coincidence with neutron capture, again we have chosen not to apply this 
veto in our full analysis. 

The shape of the background is the well known Michel spectrum, as it results from the 
decay of muons into electrons. It is identical to our /i + DAR V spectrum. However, the 
signal V e are detected via IBD which creates e + of energy 1.3 MeV lower than the original 
antineutrino, and so in fact the signal spectrum is 1.3 MeV lower than the invisible muon 
background spectrum. Nonetheless, given the energy resolutions of SK and HK, this shift 
is of limited use in distinguishing the signal and background. The normalization of the 
invisible muon background has been well-measured by SK [23] although, as these neutrinos 
are ultimately generated by cosmogenic muons, the rate may be time-dependent at the 10 - 
20% level. In summary, in combination with some accelerator down-time, not only the shape 
but also the invisible background muon rate will be known rather precisely. 

In our calculations, we have assumed that the normalization of the invisible muon back¬ 
ground and v e backgrounds are known with 7% and 10% precisions respectively, reflecting 
statistical fluctuations in the SK sample [23] and, since SK measured the total background, 
we assume an error correlation of —0.24 between the two backgrounds. 

The CC events that produce the muon also produce other particles. In fact in the case 
of A resonance CC events, which account for most of the events in the second column 
of Table 2, additional pions and generally other particles are always produced, such as 7 
[53]. These other particles in general produce some effect which is visible at the detector, 
allowing for vetoes of the background. We have studied the potential veto efficiencies. In 
our main analysis we use the selection of vetoes summarized in Table. 1. We make the crude 
approximation that all 7 events, except for those resulting from neutron capture, can be 
identified. 

For example, consider the following event. A A and an invisible /j, are created in a 
CC event. The (i decays after a few microseconds and produces an electron, leading to 
an electron-like ring. Let us call this ring number one. The A decays, for example, into a 
charged 7r. The 7 r decays quickly, yielding a /i which decays after a few microseconds yielding 
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Proposed Selection Cuts 

All events that occur during the J-PARC beam spill 
Events with multiple rings 

Events with v energies outside the 30-55 MeV window 
Events which are not fully contained 
Events within 1 ms of electron or muon rings 
Events followed by 7 emission within 1 ms 

Table 1: Proposed background selection cuts 

another e and so a second electron-like ring, let us call it ring number 2. Note that in this 
case, as the two electrons are created from the decays of distinct muons, which are created 
essentially simultaneously, the expected time difference between rings number one and two 
is roughly the muon lifetime. This background can essentially be eliminated by removing 
events with two electron-like rings. We veto events with multiple electron-like rings. 

We have folded the results of GENIE simulations of atmospheric neutrino events in water 
with the atmospheric neutrino fluxes of Ref. [27] , oscillated using the neutrino mass matrix 
parameters of Subsec. 2 . 1 , to determine the veto efficiencies. On the other hand, we use 
Super-K measurements to fix the overall normalization of the invisible muon background 
which we use in this section. GENIE yielded the number of invisible /i + and /i _ . Essentially 
all of the /i + and 80% of the /i“ come to rest and then decay in water, and so the number 
of background events is the number of n + plus 80% of the number of /i~ events. 

We End that few invisible muon events arise from NC interactions, most of which produce 
additional particles which can be used to veto them. Therefore in our main analysis we only 
consider the CC invisible muon background. In Table 2 we provide the number of invisible /i 
events caused by interactions of v and v per six years. The first row includes all events. The 
second row provides the number of events in which, according to our GENIE simulations, no 
7 is emitted as the struck nucleus de-excites. Note that 7 may also be emitted by neutron 
capture on H, however the energy of this 7 is lower than that which accompanies a nuclear 
de-excitation in most of our simulations. Furthermore we suspect that the 7 emission from 
nuclear de-excitation is generally much faster than 7 capture on H and so these may be 
distinguished. If a detector is loaded with 0.1% Gd then the vast majority of n captures 
will be on Gd. This capture results in the production of 3-5 7 which have similar energy 
distributions to de-excitation 7 , however the n capture is delayed by 10 s of /rsec and this 
delay can be used to discriminate the two kinds of 7 . 

In our analysis we use the second row, corresponding to a veto of events with additional 
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CC: E v <300 

CC: E v >300 

NC: E v <300 

NC: E u >300 

all inv. n 

343/135 

83/12 

0/0 

38/20 

no 7 

149/73 

20/2 

0/0 

17/8 

1 n 

54/99 

13/5 

0/0 

21/12 

1 n, no 7 

30/58 

5/2 

0/0 

10/5 


Tabic 2: Number of v^/v^ invisible // events per 6 years expected at SK including various 
veto conditions. The columns from left to right correspond to CC events with neutrino 
energies less than and greater than 300 MeV and NC events with neutrino energies less than 
and greater than 300 MeV. 

7 . We restrict our attention to neutrinos with E u < 300 MeV because we find that higher 
energy v essentially always create extra rings which can be used for a veto. According to the 
second row of the first column of Table 2, this leaves us with 149 invisible fjT events and 73 
invisible p + events, for a total of 222 events. Furthermore, to reduce the spallation, diffuse 
supernova backgrounds and in particular quasi-elastic interactions, we restrict our attention 
to events in which the final e energy yields a reconstructed v energy of 30-50 MeV. This 
latter condition leaves 155 of the 222 original invisible /j events, as can be seen in Fig. 3. 

In the third row no condition is placed on 7 but we consider only events which yield 
precisely 1 n. As we do not consider Gd doping, the efficiency with which such a veto may 
be implemented will be limited and so it is not considered in our analysis. Finally in the 
last row we impose both the single n and the no de-excitation 7 requirements. 

As SK had a hard trigger during the SK runs used in Ref. [23], it could not yet detect 
many of the low energy 7 ’s used in the various cuts in Table. 2. Thus it is not surprising 
that the background rate observed at SK is between the total and no 7 rates in Table. 2. 

4 Sensitivity to 5 

4.1 Optimizing the Baselines 

Where should the pDAR source be placed? As the event rate will be lower at a far detector, 
it is reasonable that the far detector be larger. This suggests that the source be to the north 
of SK. Furthermore, a maximum synergy is achieved when the baselines to SK and HK are 
as different as possible. The maximum difference is the distance between the Mozumi and 
Tochibora mines, 8 km. Therefore we will always assume that the baseline to HK is 8 km 
larger than that to SK, and will optimize the baseline to SK. 


13 












Event Type 

SK at 15 km 

HK at 23 km 

IBD Signal 5 = 0 ° 

298 

3419 

IBD Signal 5 = 90° 

325 

4549 

IBD Signal 6 = 180° 

240 

3556 

IBD Signal 5 = 270° 

214 

2426 

Invisible /i Background 

155 

3862 

Atmos v e Background 

26 

639 

p~DAR Background 

2.7 

27 

CCQE v e — 16 0 Background 

32 

332 

v f; — e~ Elastic Background 

46 

478 


Tabic 3: Total number of signal and background events expected in 6 years, implementing 
the vetoes in Table 1. 

We consider three different cases. In each case a 6 year pDAR run is combined with 12 
years of T2K (and 6 years of T2HK) in v mode only, as we find that the time in V mode 
reduces the performance of TNT2K. In the first case, only SK is considered. In the second 
only 20% of HK will be in operation for 6 years. Recall that HK consists of 10 identical 100 
kton modules, each with a 56 kton fiducial volume. Therefore, 20% of HK corresponds to 
the construction of two modules. We refer to this case as HK/5. In the third, the full HK is 
assumed to be in operation. In each case we assume that the neutrino mass hierarchy, NH 
or IH, is known but we consider both hierarchies. 

We calculate, for each baseline to SK, for 5 tme = 0°, 90°, 180° and 270°, the fit values of 
<5 fit for which the best fit of one theoretical data set to the other yields y 2 = 1, corresponding 
to an expected Ay 2 = 1 for a fit of real data. As expected this occurs for one value of <5fi t 
which is larger than <5 tm e and one which is smaller. These intervals are roughly symmetric 
about the true value of <5. In Fig. 5 we report half of the size of the interval, which is 
approximately the lcr precision which can be expected in a measurement of <5 at TNT2K 
with SK only. 

One sees that the optimal baseline depends on sin(<5). For sin(<5) = 0 it is 15-20 km, 
where the expected uncertainty on <5 is about 14° whereas for maximal CP violation it is 
20-30 km, where the expected uncertainty is about 34°. If 15 km is adopted then one Ends 
that the uncertainty lies in the range 14° — 36°, not far from the optimal for any value of 
5. However we have also found that, while 15 km yields a competitive measurement of 5, a 
longer baseline would yield a more robust breaking of the 5 —>■ 180° — 5 degeneracy. 

In Fig. 6 we present the precision with which 5 can be measured in the cases with one 
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Figure 5: The precision with which 5 can be determined using SK only as a function of the 
/iDAR baseline. The precision quoted is the average of the upper and lower uncertainties. 
Both hierarchies are considered, but it is assumed that the hierarchy is known. 
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fifth of ffK (the upper panels) and all of ffK (the lower panels) , for a given baseline from 
the /iDAR source to SK. Not surprisingly, as ffK is 8 km further, the optimal baseline to SK 
is now shorter than the SK only case. In general the most precise determination of 5 occurs 
for a baseline of roughly 15 km to SK and 23 km to HK. With one fifth of HK (all of HK) 
the precision with which 6 can be measured ranges from 9° (7°) for no CP violation to 20° 
(11°) in the case of maximal CP violation. 

This can be compared with the performance of T2HK without /iDAR and with the J- 
PARC beam running for 1.5 years in v mode and 4.5 years in v mode, with the full HK. In 
that case one expects to measure 5 [3] with a precision of 9° — 24°: With / lDAR and one 
fifth of HK, one can determine 5 more precisely than with all of HK and no pDAR. 

4.2 Measuring 5 with TNT2K 

In Fig. 7 we plot the y 2 value of the best fit of the <5fi t theoretical spectrum to the theoretical 
spectra of <5 true for SK only, at 15km and also at 23 km assuming the normal hierarchy. 
Looking at the x-axis for the value of 5 true , one observes that for a 23 km baseline maximal 
CP-violation, corresponding to 6 = 90° or 270° can be distinguished from no CP-violation, 
corresponding to 5 = 0° and indeed also to 180°, at about 4a of confidence and in fact nearly 
5 cr for 5 = 270°. With a 15 km baseline, 5 = 0° can be distinguished from S = 90° and 
270° with a bit under 4a and 6cr respectively. In both cases, 5 = 0° and 180° can only be 
distinguished at 2 — 3a of confidence. 

In Fig. 8 we consider SK at 15 km and one fifth or all of HK at 23 km. With just one 
fifth of HK, one sees that maximal CP violation, 6 = 90° (270°) can be distinguished from 
6 = 0 at more than 6a (7a). Also 5 = 0° and 6 = 180° can be distinguished with nearly 
6 a of confidence. Thus, the 5 —> 180° — 5 degeneracy, for large CP violation, is completely 
broken already with one fifth of HK. On the other hand, with the full HK, these distinctions 
can be made at more than 9a. Thus in general the full HK only serves to provide a precise 
determination of S, one fifth of HK is quite sufficient to qualitatively understand leptonic 
CP violation. 

5 Conclusions 

A 7 MW proton accelerator, with a proton energy of 600 MeV-1.5 GeV, when striking a 
medium to high Z target creates 30-50 MeV V M via p + decay at rest. We have advocated 
placing such an accelerator 15 km north of SK. The oscillated V e can be detected by SK via 
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NH T2Kv+T2HKv+(o.DAR (5%) IH 




Figure 6: The precision with which <5 can be determined using SK and one fifth of ffK (top) 
or all of ffK (bottom) as a function of the //DAR baseline to SK, as in Fig. 5. The baseline 
to ffK is 8 km greater. 
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X 2 for SK + ^DAR [15 km] 



X 2 for SK + uDAR [23 km] 



Figure 7: y 2 value of each trial value of <5 vs each true value assuming the normal hierarchy, 

using SK only at 15 km (top) and 23 km (bottom) with 6 years of /xDAR and 12 years of 

T2K operating in v mode. From the .x-axis one observes that maximal and null CP-violation 

can be distinguished at about 4 — 5cr while, with null CP-violation, the sign of cos(<5) can be 

determined with 2 — 3a of confidence. Null CP violation can be excluded at more than 3a 
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of confidence for more than half of the values of 5. 




X 2 for SK + HK/5 + |xDAR [15 and 23 km] 



X 2 for SK + HK + |xDAR [15 and 23 km] 



Figure 8: As in Fig. 7 but now including HK (bottom) and also one fifth of HK (top). With 
one fifth of HK, maximal and null CP-violation can be distinguished at more than 7rx while, 
with null CP-violation, the sign of cos(<5) can be determined with nearly 6cr of confidence. 
Null CP violation can be excluded at more than 5rr of confidence for more than half of the 
values of 5. 
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IBD. The V spectrum is known quite precisely as is the IBD cross-section, thus systematic 
errors are small in the determination of the CP-violating phase 5. Together with 12 years of 
T2K, this setup allows 5 to be measured with a precision of 14° — 36° in 6 years. Maximal and 
null CP violation can be distinguished at about 4 — 5cr and 5 = 0° can be distinguished from 
5 = 180° with 2 — 3<v of confidence, a very difficult task for conventional beam experiments. 

If just one fifth of HK is built, corresponding to two of the ten planned modules, then 
at the preferred Tochibora mine it will be 23 km from the yuDAR source. This will allow an 
excellent determination of S, with a precision of 9° — 20° in 6 years. Maximal and minimal 
CP violation can be distinguished with 6 — 8a of confidence and the sign of cos(J) at well 
over 5 a. 

The required accelerator is not far beyond the state of the art, and as the neutrinos 
are created from decay at rest the only requirement is that the tt and /i stop in the target. 
This means that many other physics programs can be done simultaneously with the same 
accelerator, for example it can run an accelerator driven subcritical nuclear reactor. Such 
accelerators in fact have just the specifications that we require. 
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